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Abstract 

We analyze here the minimal conditions for directional motion (net flow in 
phase space) of a molecular motor placed on a mirror-symmetric environment 
and driven by a center-symmetric and time-periodic force field. The com- 
plete characterization of the deterministic limit of the dissipative dynamics of 
several realizations of this minimal model, reveals a complex structure in the 
phase diagram in parameter space, with intertwined regions of pinning (closed 
orbits) and directional motion. This demonstrates that the mirror-symmetry 
breaking which is needed for directional motion to occur, can operate through 
an internal degree of freedom coupled to the translational one. 
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The possibility of motion rectification at brownian scale from non equilibrium fluctuations 
is an interesting question already posed by Smoluchowski |I| and Feynman 0. It has 
been shown that directional motion (appearance of net flow in phase space) of brownian 
particles, without any macroscopic gradient or force, can be achieved provided that the 
potential exhibits spatial asymmetry and detailed balance is broken. This latter condition 
assures that the system is out of equilibrium. Systems with a periodic potential profile but 
spatial asymmetry, called ratchets, have been addressed as systems in which non equilibrium 
fluctuations can induce directional motion ||. They have attracted much attention on the 
basis that they can help to understand the physics of molecular motors 0, along with the 
possibilities they open for using these ideas in superconductors 0, Josephson junctions ||, 
quantum dots |7j and in the promising world of nanotechnology [|J. 

Hereon, we will focus on "motor" systems such that their internal degrees of freedom are 
essential for net directional motion to occur. The motivation for this problem came initially 
from the (bio)molecular motors field when it was found || that kinesin direction of motion 
along microtubules could be reversed by modifying the architecture of a small domain of 
the protein called "neck region". These discoveries could suggest that the mirror symmetry- 
breaking mechanism responsible for the directional motion performed by these proteins could 
lie in their own structure rather than in their environment. In this paper we will consider the 
type of systems with a mirror-symmetric environment for its traslational degree of freedom 
(u tr ). The symmetry-breaking mechanism acts through an internal degree of freedom (u mt ) 
coupled to it, that is, we will consider a dimer model (two degrees of freedom). We are 
interested in the minimal conditions for the operation of this system as a motor (being 
capable of moving against an applied field). Previous works have also analysed directional 
motion in dimer models |T0[ but the potential for u tr considered was ratchet ab initio. On the 
contrary, the dimer model we present in this paper is inmersed in a symmetric environment. 
This is also the case for the system recently studied by |ITJ|. There, internal degrees of 



freedom (more than two) experience a flashing interaction potential, while ours is a rocked 
system. After a complete characterization of the deterministic limit of the dynamics, which 



reveals the basic nonlinear mechanics of directional motion, we will briefly discuss the gross 
features of the Langevin dynamics of this model and its utility as a help to design new 
technologies at micro and nanoscales. 

As a working mechanical visualization of this model, one can use two overdamped and 
coupled brownian particles, with positions u\ and u 2 (so u tr = 1/2(m 1 +m 2 ) and u int = U2—U1) 
moving in a periodic, symmetric potential V(u + 1) = V(u) = V(— u) and being driven by 



center-symmetric periodic forces Fi = —Fi(t + T/2)i = 1,2 of period T = 2tt/w \ 12j . In the 
deterministic limit, the equations of motion read 

ii x = -V'(u 1 )-d 1 W(u l ,u 2 )+F 1 (t) (1) 

u 2 = -V\u 2 ) - d 2 W{ Ul , u 2 ) + F 2 {t) (2) 

We impose on W(u\,u 2 ) the general condition of being a function of the relative distance 
u mt = u 2 — Mi, so the partial derivatives with respect to u\ and u 2 verify diW{ui 1 u 2 ) = 
— d 2 W(ui, u 2 ). We will consider the cases in which W is a convex and a non-convex function 
of u mt . 

Equations (l)-(2) remain invariant under the symmetry transformations (ui,u 2 ,t) — > 
(-tta, -ui, t+T/2) provided that F x (t) = F 2 (t) and (u u u 2 ) -> (-u 2 , -Ui) ifF^t) = -F 2 {t). 
With this proviso one can easily show that any averaged "directional motion" in phase space 
is necessarily (and straightforwardly) null : if there exists a solution of (l)-(2) with nonzero 
velocity, by symmetry we can find another solution with the same velocity but opposite 
sign, so no net motion can be observed when averaging over the phase space. Directional 
motion in this strong sense can only occur if the inequality Fi(t) 7^ ±F 2 (t) holds, regardless 
the specific form (convex or non-convex) of the interaction potential W. Taking Fi(t) from 
the class of functions F(t) = F ac sin(co>t + 2n5) (the simplest periodic center-symmetric 
function F(t) = —F(t + T/2)), there are two ways of breaking the symmetry of this system 
(F 1 (t) ^ ±F 2 (t)): 

• i) max t Fi(t) 7^ max t F 2 (t) 



• ii) 5\ 7^ 5 2 mod | 

That is, applying forces on each particle with (i) different amplitude or (ii) different phase. 
We will focus on the limit cases for each situation: we will consider F^) (t) = and Fj® (t) 7^ 
as well as 5\ = and 5 2 7^ 0. Any other situation verifying the condition Fi(t) 7^ ±F 2 (t) 
is just a combination of these two possibilities. 

In particular, most numerics have been performed with the simplest choices 

V ( u ) = 77T^cos(27ra) W(u u u 2 ) = ^[(u 2 - u x ) ~h] 2 (3) 

The convexity of W and the dissipative character of the dynamics have the important con- 
sequence that the partial order among initial conditions is preserved (monotonicity property 
or "no-passing rule" [|13|]). The monotonicity property says that if at a time t = t two 



initial conditions u,v verify that - u 1 (t ) < ^1(^0) an d ^2(^0) < ^2(^0)) then for any time 
t > to the inequalities u\(t) < vi(t) and u 2 (t) < v 2 (t) hold. This property leaves small room 
for deterministic complexities of this non-integrable dynamics. For example, the asymp- 
totic mean velocity of all trajectories in the phase space is unique, and vibrating pinned 
solutions (closed orbits) cannot coexist with mobile ones. The choice of W(u\, u 2 ) im- 
plies another symmetry relation: the system remains invariant under the transformation 
(ui,u 2 , Zo) — > (—ui, —u 2 , 1 — lo). That is, if we have a mobile solution of (l)-(2) for a value 
Iq we have another mobile solution with opposite velocity for 1 — Z - For fixed values of 
(Q,K,l ) we have approximated numerically to optimal accuracy the function J(F aci uj) 
where J is defined as the asymptotic flow in phase space J = (u tr ). 

In figures l.a) and l.b) the J(F ac ) profiles at typically low (u = In x 0.01) and high 
frequencies (u = 2tt x 0.1) for Fi(t) = and F 2 (t) = F ac sin(ujt) (case i) — different 
amplitudes — ) are shown for the values (Q,K,l ) = (2, 1, 1/4). As can be seen from figure 
l.a, the complex stairlike structure of the relation J(F ac ) for low frequencies is qualitatively 
indistinguishable from those exhibited by single particle rocking ratchets (asymmetric po- 
tential). However, when the frequency is increased, the direction of motion for fixed Iq can 
be either positive or negative (figure l.b), something which cannot occur in deterministic 



single particle rocking ratchets. The primary structure of mobility bands in F ac separated by 
vibrating pinned solutions is well understood in terms of saddle-node bifurcations signaling 
the onset of global flow in phase space. Figure l.c shows the bifurcation diagram of pinned 
(vibrating) solutions. <p represents the mean value over a period of the external force of 
the traslational variable, = (u tr )T- Thick lines represent stable orbits whereas thin ones 
represent unstable ones. This bifurcation diagram neatly determines the intervals of motion 
in F ac as well as the local sign of J in both interval's edges. The numerical inspection of the 
scaling of intermittencies density is the one of a type I intermittency scenario (associated 
here to quasiperiodicity and frequency locking), giving the staircase aspect of J(F ac ) close 
to the mobility onset. 

In figure 2. c an example of jT(F ac ) profile for F x (t) = F ac sm(ujt), F 2 (t) = F ac sm(ut+27i5) 
(case ii) — different phase — ) with (u,5) = (2tt x 0.05,0.35), (Q,K,l = (2, 1/2, 1/4) shows 
clearly the complex alternation of positive and negative flow in the phase space. Again 
the mobility bands in F ac when (a>, 5) remain fixed, appear after the collision of stable and 
unstable orbits by means of a saddle-node bifurcation. In figure 2. a) we have plotted for 
uj = 2tt x 0.05 the width in F ac of the mobility bands at different values of 5. Figure 2.b) 
shows the values of F ac at which the first and second bifurcations occur, that is, the smaller 
value of F ac at which the stable and unstable orbits collide (depinning transition) and the 
nearest value at which they emerge again (pinning transition) what determines the width of 
the first mobility band. As usual, very simple dynamics reveal astonishingly complex phase 
diagrams whenever time and length scale competition plays a role. 

We want to stress that the rectification mechanisms in the models presented here are 
completely different from other rocking ratchet systems with internal degrees of freedom 
(d.o.f) |TIJ, where the symmetry is broken ab initio by the asymmetric potential V(u). The 
only system with a symmetric environment for the traslational d.o.f. in which symmetry 
breaking comes through the internal d.o.f. is the flashing system (as oposed to rocking) 
studied by Porto et al. [TTJ|. We also remark that the rectification mechanisms for the cases 



numerically studied in the preceeding paragraphs show some differences: in case (i) we 



observe directional motion even at the adiabatic limit (slow varying forces) whereas in case 
(ii) directional motion occurs only at finite frequencies. When applying forces with different 
amplitudes, the ratchet effect lies in the fact that one particle acts as a cargo, so it is easier 
to move the system in the " driver" particle direction than in the other. In figure 3 it can be 
clearly seen that the depinning force is smaller in the driver's direction than in the other. 
When applying different phases, the ratchet effect is more subtle: the combined effect of 
phase value and strength F ac determines which particle plays the role of cargo and which 
one the role of driver. The absence of the adiabatic limit for this case (ii) precludes the use 
of time independent schemes for understanding rectification in an intuitive way. In figure 
4 two trajectories corresponding to both limit situations (cases i), ii) above) analyzed are 
drawn. Although our system is a rocking ratchet, there is an alternative view in which, 
looking at the equations of motion for the traslational variable u tr 



one can see for both situations, case i) and ii), that the time dependence (periodic or 
quasiperiodic) of u mt allows the interpretation of the first term in rhs as a flashing potential 
for u tr . The directional motion can thus be seen as the result of the adequate synchrony 
between external force on u tr and the periodic flashing potential. 

When the convexity condition on the interaction potential W(ui,u 2 ) is removed and 
non-convex interaction (such as double- well or Lennard- Jones) potentials are considered, the 
monotonicity property (which severely restricts the complexity of the dynamics) is lost, and 
generically the space is partitioned in basins of attraction corresponding to attractors with 
different asymptotic velocities. The description of the dynamics becomes thus more complex. 
The phase portrait shows, for certain regions of parameter values chaotic attractors and there 
appear new bifurcations (other than the observed for the convex case: pitchfork, generic 
saddle-node and stability interchange) as the parameters change. Anyway the existence of 
net flow (ie. non zero phase space average of asymptotic velocities) in phase space is a 
generic property in wide regions of parameter space, as in the convex case. 
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We are dealing with rocking-like systems where the rectification mechanism is com- 
pletely deterministic, so when introducing noise £i(t) and £ 2 (£) m equations (l)-(2) (white 
gaussian noise with correlation function (£i(t)£j(s)) = D6ij6(t — s), i,j = 1,2, being 
D = ksT the diffusion coefficient) the mobility bands widen (fluctuation-induced depin- 
ning) and J decreases with increasing noise strength (within the deterministic mobility 



bands). Both phenomena are easily explained in the limit of small noise [JT4] . For high D 
the phenomenon of current reversal arises too and when noise strength is high enough 
diffusive motion dominates and no directional motion can be observed. 

The proteins of the kinesin superfamily are just composed of two globular "heads" 
(catalitic domains that move over the microtubule using the energy delivered in ATP- 
hydrolysis). As mentioned before, recent experimental work suggests that the mechanism for 
directionality may rest on the characteristics of the flexible structure joining both domains. 
ATP-hydrolysis induces conformational changes in these proteins, that is, acts through their 
internal degrees of freedom what may provide the symmetry-breaking needed for the di- 
rectional motion along the microtubule. The detailed modelling of these conformational 
changes of the kinesin is beyond the scope of this work; nevertheless the generic model dis- 
cussed in this paper could help to understand in the simplest mechanical terms the role of 
the internal degrees of freedom in molecular motors. 

Rigorous results on general simple models like this, could serve to guide or inspire new 
technological applications specially in the new field of nanotechnology. A simple nano-scale 
realization of the type of motors considered here would consist of clusters of entities (particles 
or macromolecular agregates) with different electrophoretic mobilities joined with the aid 
of "flexible" molecules or polymers. In order to reproduce the symmetric environment it 
will suffice to place this engine in a row of symmetric obstacles (electrodes for instance). 
When applying an ac electric field, because of the difference in the electrophoretic mobilities, 
different periodic forces will act on each cluster, making it possible to observe directional 
motion. 

In summary, the analysis of these minimal models, convincingly demonstrate that the mirror 



symmetry breaking needed for directional motion to occur, can act through an internal 
degree of freedom eventhough the overall position of the system experiences a symmetric 
environment. 
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FIG. 1. a) and b) Flow as a function of the amplitude of the external force F ac at two different 
frequency values, c) Bifurcation diagram using the mean value over a period of u tr 4> = {u tr )T as 
relevant magnitude for u = 2w x 0.1. Intervals between vertical lines correspond to mobility bands 



in F ac . 
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FIG. 2. a) Mobility bands in F ac at u; = 2ir x 0.05 as function of the phase difference 5. b) 
Lower and upper values of F ac for the first mobility band as a function of 5. c) Dependence of the 
velocity with F ac for (u,S) = (2vr x 0.05,0.35). 
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FIG. 3. Flux J in the presence of a constant force: it can be clearly seen that the depinning 
force is smaller in the "driver" particle direction (positive one). 
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FIG. 4. Deterministic trajectories (u\(t) dashed line and U2(t) solid line) for the two models 
discussed, a) corresponds to different amplitudes (case i) (uo,F ac ) = (2ir x 0.01,0.35) and b to 
different phases (case ii) (u,F ac ,S) = (2ir x 0.05,1.284,0.35). 
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